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Abstract The static behavior of orientational glasses is discussed in terms of 
a replica theory based on the infinite range random bond-random held model. 
A general version of the model applicable to dipolar and quadrupolar glasses 
is presented using a symmetry-adapted representation for the order parameter 
fields. Numerical results for the < 100 > quadrupolar glass are obtained and 
compared with the Ising and < 111 > orientational glass models. 



1) INTRODUCTION 
Q 

qq Orientational glasses 1 (OG's) have been attracting considerable experimental and 

t-h theoretical interest since they are believed to provide a conceptual link between spin 

glasses and 'canonical' real glasses. A glass-like transition in these systems is due to 
^ the freezing-in of molecular orientations into a disordered configuration. Physically, 

— molecular orientation is associated with either an electric dipole or an electric or 

elastic quadrupole moment, thus one may also speak about dipolar and quadrupolar 
^ glasses, respectively. Typical examples of dipolar glasses are the 'proton' glass 2 

ON Rbi_ x (NH4) x H 2 P04 or RADP and related compounds. Among quadrupolar glasses, 

mixed alkali cyanide crystals KB ri _ x (CN) x 3 ' 4 , N ai _ x K x CN 5 ' 6 ' 7 , Na(Cl)i_ x (CN) x 6 ' 8 
etc. have been the subject of several recent investigations. While dipolar and 
^ quadrupolar glasses share many analogous features, there exist some differences due 

to their symmetry properties. For example, dipolar glasses possess a held inversion 
O symmetry, which is absent in the quadrupolar case 9 . 

In the present paper we discuss a simple theoretical model believed to be 
capable of describing the mechanism of the OG transition in dipolar and quadrupolar 
glasses. In particular, we will focus on quadrupolar glasses of the mixed cyanide 
type in which the rotation of the CN~ ions is strongly hindered by the crystal 
anisotropy. This then restricts molecular orientations to a set of discrete directions, 
which lie along the < 100 >, < 111 >, or < 110 > crystal axes in an effectively 
cubic crystal. The corresponding discrete-state model has recently been derived by 
semimicroscopic arguments by Vollmayr et al. 10 and was shown to be equivalent to 
the 5-state Potts model where 5=3,4,6, respectively, for the above three cases. This 
can formally be regarded as a generalization of dipolar glass models, where an electric 
dipole is typically assumed to have just two discrete orientations corresponding to 
two equilibrium states in a bistable potential. Such systems can be described by 
an Ising pseudospin model 11 or equivalently the 2-state Potts model. It should be 
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noted that the s = 2 pseudospin model is generally not applicable to quadrupolar 
glasses except in some special situations, for example, when molecular rotation is 
restricted to two discrete perpendicular orientations in a plane and the external 
fields are applied within the same plane. 



SYMMETRY-ADAPTED ORDER PARAMETER FIELDS 

In a discrete-state model the orientational degrees of the z-th molecule are repre- 
sented by a set of occupation numbers Ni p = 0,1, where p = labels the 
orientations, which obviously satisfy the sum rule 

E^ P = i, (i) 

P =i 

with s=2 for dipolar and s = 3,4,6 for quadrupolar glasses. 

In constructing a realistic model of an OG one should account for some charac- 
teristic symmetry features not found in magnetic spin glasses. Namely, experiments 
show that the OG order parameter is non-zero in the high temperature phase, sug- 
gesting the existence of random fields 12 acting on the orientational degrees of free- 
dom in addition to random pairwise interactions of the spin-glass type. This is the 
essence of the random bond-random field (RBRF) model of dipolar glasses 11 , which 
can easily be generalized to the case of quadrupolar glasses by including the appro- 
priate random strain terms into the quadrupolar glass Hamiltonian 13 . In the mixed 
cyanide systems, random strains can be uniquely represented by their irreducible 
components with A\ gi T 2g} and E g symmetries. Similarly, the orientational degrees 
of freedom can be expressed as linear combinations of occupation numbers Ni p , 
which will transform according to the same irreducible representations. By analogy, 
in continuous models of OG's the orientational degrees of freedom are described in 
terms of symmetry-adapted spherical harmonics i^(t?, <£>),(// = 1, 5), where the 
angles i?, <p specify the orientation of the molecular axis 12 ' 14 . Thus we will introduce 
new symmetry-adapted order parameter fields 13 

s 

Zifi = A^pNip , (2) 
P =i 

where // now labels the appropriate irreducible represenation. The coefficients A w 
are simply determined by group theory arguments. Dropping the site indices i and 
labeling the orientations p in the same way as done earlier 10 ' 13 , we find for the three 
representative cases of CN~ orientations: 

< 100 > : Z 1 = y/3/2 (TVi - N 2 ) ; Z 2 = ^1/2 (2iV 3 - N t - N 2 ) . (3) 



< 111 > 



z 2 

Zs 



= N 1 + N 2 -N 3 -N 4 ; 
= N 2 + N 3 -N 1 -N 4 ; 
= N 3 + N 1 -N 2 -N 4 . 



(4) 



< 110 > : Z x = (TVi - iV 4 ); Z 2 = ^3 (iV 2 - iV 5 ); Z 3 = y/Z (N 3 - 7V 6 ); 
Z 4 = ^Jz/2(N 2 + N 5 -N 3 - N 6 ) ; 

Z 5 = yJT/2(2N 1 +2N 4 -N 2 -N 3 -N 5 -N 6 ) . (5) 

It should be noted that the number of non-trivial fields is in each case equal 
to 5 — 1, and the remaining field is given by Z s = 1 in view of the closure relation 
(1). The symmetry- adapted fields in Eq. (4) and the first three fields in Eq. (5) 
transform according to the irreducible representation T 2g} whereas the remaining 
fields belong to the E g symmetry. 

The case s = 2 in Eq. (2) corresponds to an Ising pseudospin S = Z\ = 
Ni — N 2 = ±1, where of course N x + N 2 = 1. The algebra of the Z^-fields for 
s > 2 can simply be derived from the occupation numbers Ni p ; it will in general 
differ from the Ising case thus reflecting the distinct symmetry character of the 
quadrupolar glass models. 

Besides reducing the number of relevant fields from s to s — 1 the symmetry- 
adapted representation leads to a simple parametrization of the quadrupolar glass 
order parameter tensor in terms of its symmetry components. In the following we 
will formally derive the OG free energy functional for s = 2, 3, 4, 6 and then discuss 
in detail the < 100 > RBRF model and compare the results with the < 111 > and 
Ising cases. 

RANDOM BOND RANDOM FIELD MODEL 

Assuming isotropic random interactions Jf, 10 ' 13 and using the symmetry-adapted 
representation Z 8|U the RBRF model of OG's can be written in the general form 

s — 1 s — 1 

'H = — 2 X! X] Zif+Z^ — ^2 ^2 hipZi^ . (6) 

ij /i=l i /i=l 

Here /j^ are random fields representing local random electric fields or random strains 
in dipolar and quadrupolar glasses, respectively. One may also add a homogeneous 
external field of symmetry // to /j^. As usual, we assume that both random 
bonds and random fields are Gaussian and uncorrelated, and can be characterized 
by their first and second cumulant averages [</ij] = Jo/N, [{Jij) 2 ] C av = J 2 /N and 
[hi»] av = 0, [hi^hjv]^ = A^SijS^, respectively. 

The model (6) with Z 8|U from the set (3) is believed to be applicable to mixed 
cyanide systems such as Nai_ x K x CN 15 or KBri_ x (CN) x 16 in the concentration range 
where the most probable orientations of the CN~ axes are along < 100 >. The case 
< 111 > has already been discussed earlier 13 . 

To obtain the average free energy we now apply the usual replica method 17 by 
adding a replica index a = 1, . . . , n to all fields Z" in the Hamiltonian and writing 

(7) 



T = -/T 1 lim 



d 



o dn 



Tr n exp 1-/3 ^ H c 



a = l 
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After performing the random averages and using the standard manipulations from 
the theory of spin glasses 17 we obtain for the average free energy per CN~ molecule 



1 - 9 l ~^Y,(P a J-\ 

+ In Tr n exp 



p4 8) £ + ft 2 j 2 E E (C + ^)z a X 

where the sum Y^( a p) runs over distinct pairs of replicas with a ^ f3. The parameter 
5 refers to a specific model, i.e., s = 2 for Ising and s = 3,4,6 for the < 100 >, 

< 111 >, and < 110 > RBRF models, respectively, and the effective average inter- 
action constant is given by 

4 S) = Jo + \{s - 2)/3J 2 . (9) 

The equilibrium values of the order parameters P^ and are obtained from the 
saddle-point equations df /dP^ = df /dq^ = 0, yielding 

P? = (ZZ), C = < z l z »)- (io) 

REPLICA-SYMMETRIC SOLUTION: MODEL < 100 > 

Next we look for a replica-symmetric solution for the order parameters (10) of model 

< 100 > (s=3) by writing = Pi, P 2 ° = P 2 , and q^ = q^ q 22 ° = q 2 , ql 2 = q 2 i = qr 
for all a ^ f3. After linearizing the exponential terms in Eq. (8) we find 



/ = -'/J' " I) 2 + (?2 " l) 2 + 2q T ] + i/Jjf (P 2 + Pi 

/~\~ OO f' ~\~ oo 

dx\ J dx 2 exp ( — x 2 /2 — x 2 /2j In g(H\, H 2 ) , (11) 
-oo J — oo 



27T 

using the notation 

g(H 1 ,H 2 ) = 2cosh(y3/2/3#i)exp + exp (2V2f3H 2 ) (12) 

with 

H x = J( qi + A - q T ) 1/2 Xl + 4 3) P 1 + (V2/2)/3 J 2 q T ; 
H 2 = J(q 2 + A - q T ) 1/2 x 2 + 4 3) P 2 + {V2/8)pj 2 { qi - q 2 ) . (13) 

The conditions df/dP\ = . . . = df/dqx = then lead to a set of coupled equations 
for the order parameters, which can be written in compact form by introducing the 
local strain polarizations pi and p 2} 

p^dlng/diPH,), (// = 1,2) , (14) 




Figure 1: Temperature dependence of the quadrupolar glass order parameter for the 
< 100 > RBRF model, for several values of the random-held strength A. 

namely, 

p n = M12 > % = [pI] 12 > ?t = [pm] 12 , (15) 

where [. . .] 12 means a Gaussian average over the variables x\ and x 2 . 

For A < 1 these equations can be linearized with respect to P^ and and it 
follows that for T > T n , where 

1/2- 



2 L 



J + J ' + 2J' 



(16) 



one has Pi = P 2 = 0. The phase diagram is similar to the one obtained for a 
continuous isotropic model of the quadrupolar glass 18 . In particular, for J < 
and | J 1 3> J there will be no long-range order down to the lowest temperatures. 
Restricting ourselves to the case P^ = with A ^ 0, where similar qualitative 
arguments apply, we find by numerically solving the remaining equations that 

q T = , qi = q 2 = q • (17) 

This corresponds to an isotropic quadrupolar glass phase with a single order param- 
eter q = q(T). Numerical solution for q(T) is shown in Fig. 1 for several values 
of the random-held strength A. It is interesting to compare this result with the 
< 111 > 13 and Ising 11 cases for the same value of A (see Fig. 2). Notice that the 
differences between the three cases appear to be quite small. 

For T < T c one expects in general an anisotropic quadrupolar glass phase 
with qi ^ q 2} qj ^ in the presence of long-range order with Pi ^ 0. However, 
it should be noted that the replica symmetric solution will become unstable at low 
temperatures and one should therefore hrst investigate the limits of its stability. 
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Figure 2: Comparison between the < 110 >, < 111 >, and Ising RBRF models for 
A = 0.01 and 0.2. 

STABILITY CRITERION AND THE GLASS TRANSITION 

In analogy to spin 19 and dipolar 11 glasses, and to quadrupolar glasses with discrete 13 
or continuous 20 orientational degrees of freedom one can investigate the stability of 
the isotropic OG phase by considering small Gaussian fluctuations p^ of the order 
parameter around the solution (17). This procedure is expected to be applicable 
even in the case when the transition is weakly hrst order. Thus we will write 

C = (? + pf)s^ (is) 

and look for the deviation of the free energy (8) from the replica symmetric value 
(11). For T > T c we have 

= \P 2 J 2 Jim ^~ J2 E E {<V<W)(7S) 
z u an ^ ( ^ ^ 

-/3 2 J 2 plZlZZZl) - {Z«X){Z2Zl)\ } pWpM . (19) 

For the < 100 > model the hrst average in the square brackets is, for example, 

{Z a XZ2Zl) = [ P l P l] ^ , (a + fl ± 7 + 6) , (20) 
using the notation of Eq. (15). 




The stability criterion requires the eigenvalues of the matrix in the curly brack- 
ets in the limit n — > to be all positive 19 . For the < 100 > case we find two relevant 
eigenvalues, namely, 

A ± = i-^{i-|«+[rfL-^[rfL±(5«+[rirfL)} • < 21 > 

Obviously A + < A_, hence the transition temperature Tj is determined by the 
condition A + = 0. 

The corresponding eigenvalue for the < 111 > model is given by 

A+ = 1 - /3V 2 {l + [pt] i23 + 4 [ PlP2P3 ] 123 + 2 , (22) 

where 13 

tanh(/3ff At ) - jX^ tanh(/3^) 

^ = MiMW ' (/i ' z/ = 1 ' 2 ' 3) ' (23) 

with Hp = J(q + A) 1 / 2 ^^, and the averages [. . .] 12 3 are now over the three Gaussian 
variables Xi } x 2} x 3 . 

In Fig. 3 we plot Tj as a function of the random-held variance A for the 
< 100 > (s = 3) and < 111 > (s = 4) models and compare it with the coresponding 
result for the Ising 11 (s = 2) case. It should be noted that for a given value of A 
the Ising RBRF model has the lowest and the < 111 > model the highest transition 
temperature, indicating that stability of the replica-symmetric solution q(T) might 
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be related to the number of equilibrium orientations s, i.e., q(T) appears to be more 
stable for models with smaller values of parameter s. 

For T < Tj replica symmetry is broken and the system is in a non-ergodic 
quadrupolar glass phase. The precise form of the Parisi order parameter function and 
its dependence on A for s = 3, 4, 6 is not yet known, however, it will be presumably 
analogous to the one of the corresponding Potts glass 9 . 

DISCUSSION 

The OG order parameter q(T) can be determined experimentally by NMR 6 ' 21 and 
by x-ray diffraction 22 or similar techniques. In the deuteron glass DRADP the 
temperature dependence of q(T) obtained from the second moment of the 87 Rb 
NMR spectrum 21 was found to be in good agreement with the predictions of the 
RBRF Ising model 11 . A similar analysis was carried out for quadrupolar glasses 6 
NaCli_ x (CN) x and Na x Ki_ x CN, where q(T) was determined from the NMR spec- 
trum of 23 Na. According to Fig. 2 the difference between the results for the < 100 > 
and Ising RBRF model is rather small and thus not easily detectable in experiments 
of the above type. A more precise evaluation of q(T) for the above systems could be 
made by means of the 14 N NMR spectrum which probes locally the orientation of 
each individual CN~ ion, whereas the 23 Na line shape contains contributions from 
a cluster of neighboring CN~ ions. 

Another test of theoretical predictions can be made by measuring the tempera- 
ture at which the splitting between the static field-cooled and zero-field cooled linear 
susceptibilities %fc and Xzfc, respectively, occurs. This was indeed performed in the 
quadrupolar glass 3 KBri_ x (CN) x and in the dipolar glass 23 DRADP. By fitting the 
observed temperature dependences of %fc and Xzfc to the corresponding theoreti- 
cal expressions one can determine the model parameters and hence the value of the 
freezing temperature Tf. The problem with this method is, however, that the value 
of the static susceptibility depends on the experimental time scale. Thus a static 
value of Tj is not directly accessible in experiments carried out on a finite time scale 
and could only be obtained by a suitable extrapolation to infinite observation times. 
Dielectric experiments show that dipolar glasses are strongly polydispersive with 
maximum relaxation time T max increasing dramatically on approaching the glass 
transition. It has recently been shown that an analysis of the dielectric spectrum 
of the dipolar glass DRADP can be greatly facilitated by introducing a so-called 
temperature-frequency plot 24 , which reveals that in the glassy regime T max diverges 
according to the Vogel-Fulcher law 

W = T exp[*7/(r-To)], (24) 

where r and U are constants and the Vogel-Fulcher temperature T marks the onset 
of the glass transition and can thus be related to the freezing temperature Tj. In 
principle, a similar analysis would be applicable to other OG's including the case of 
quadrupolar glasses. 

Finally, an interesting phenomenon not yet fully understood is the behavior 
of the nonlinear susceptibility X3 of OG's, which was found in both the dipolar 



glass 25 Ki_ x Na x Ta03 and quadrupolar glass 4 KBri_ x (CN) x to obey a scaling law 
X3 ~ (T — Tg) -1 . Here 7 1.7 for the dipolar glass with x = 0.2, and 7 0.95 and 
7 1.7 for the quadrupolar glass with x = 0.41 and x = 0.53, respectively, and T g is 
expected to be a measure of Tf. It has been argued using a dynamic version of the 
RBRF Ising model that in dipolar glasses 26 the dynamic OG susceptibility Xz{ u ) 
approaches a scaling form of the above type in the limit u — > with a mean-field 
exponent 7 = 1. A generalization of this result to the quadrupolar glass case is still 
missing. As shown by Haas et al. 27 for the 3-state Potts glass in uniform and random 
fields the second-order static nonlinear susceptibility % 2 is non-zero due to a lack 
of spin inversion symmetry and diverges to —00 at the zero temperature transition, 
while X3 diverges to +00 as in spin glasses. One should try to formulate a dynamic 
theory of nonlinear OG susceptibilities for the present symmetry-adapted RBRF 
model in cases s = 3, 4, 6 and compare its predictions with the experimental results. 

CONCLUSIONS 

The RBRF of OG's has been presented in a general form which is expected to 
be applicable to both dipolar and quadrupolar glasses. The essential step is the 
use of symmetry-adapted OG order parameter fields, which transform according to 
the irreducible representations of the underlying symmetry group of the system. 
This was done explicitly for the Ising case (s = 2) of a dipolar glass and for cubic 
quadrupolar glasses of mixed alkali cyanide type with CN~ equilibrium orientations 
along the < 100 >, < 111 >, and < 110 > axes (s = 3,4,6), where s is the number 
of discrete orientations. The advantage of the symmetry-adapted representation for 
the s = 3,4,6 models is that the OG order parameter tensor can be parametrized 
in terms of its symmetry components. Explicit numerical solutions for the replica- 
symmetric OG order parameter q(T) have been obtained for the < 100 > and 
< 111 > 13 models in a broad temperature range, suggesting that these models could 
readily be used in analyzing the experimental data for real quadrupolar glasses with 
appropriate symmetries. The freezing temperature Tj separating the ergodic high- 
temperature OG phase from the nonergodic low-temperature phase can be calculated 
numerically for each model as a function of the random-field strength A. It turns 
out that at fixed A replica-symmetric solution is less stable for larger values of s, 
thus resulting in a higher freezing temperature Tj of the corresponding models. 
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